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At present cosmology provides the nominally strongest constraint on the masses of standard 
model neutrinos. However, this constraint is extremely dependent on the nature of the dark energy 
component of the Universe. When the dark energy equation of state parameter is taken as a free 
(but constant) parameter, the neutrino mass bound is m v 5= 1-48 eV (95% C.L.), compared with 

m v < 0.65 eV (95% C.L.) in the standard model where the dark energy is in the form of a 
cosmological constant. This has important consequences for future experiments aimed at the direct 
measurement of neutrino masses. We also discuss prospects for future cosmological measurements 
of neutrino masses. 



In the past few years a new standard model of cos- 
mology has been established in which most of the energy 
density of the Universe is made up of a component with 
negative pressure, generically referred to as dark energy. 
The simplest form of dark energy is the cosmological con- 
stant, A, which obeys Pa = — PA- This model provides 
an amazingly good fit to all observational data with rel- 
atively few free parameters and has allowed for stringent 
constraints on the basic cosmological parameters. 

The precision of the data is now at a level where ob- 
servations of the cosmic microwave background (CMB), 
the large scale structure (LSS) of galaxies, and type la 
supernovae can be used to probe important aspects of 
particle physics such as neutrino properties. Conversely, 
cosmology is now also at a level where unknowns from 
the particle physics side can significantly bias estimates 
of cosmological parameters. 

The combination of all currently available data from 
neutrino oscillation experiments suggests two important 
mass differences in the neutrino mass hierarchy. The so- 
lar mass difference of Am 2 2 ~ 7 x 1CU 5 eV 2 and the 
atmospheric mass difference Amf.3 — 2.6 x 10 -3 eV 2 0. 
In the simplest case where neutrino masses are hierar- 
chical these results suggest that mi ~ 0, mi ~ Am so iar, 
and m.3 ~ Am a tmosphoric- If the hierarchy is inverted one 
instead finds 7713 ~ 0, rri2 ~ Am a tmospheric, and mi ~ 
Am a tmosphoric- However, it is also possible that neutrino 
masses are degenerate, mi ~ m-2 ~ m.3 S> Am a t m0 spheric- 
Since oscillation probabilities depend only on squared 
mass differences, Am 2 , such experiments have no sen- 
sitivity to the absolute value of neutrino masses, and if 
the masses are degenerate oscillation experiments are not 
useful for determining the absolute mass scale. 

Instead, it is better to rely on kinematical probes of 
the neutrino mass. Using observations of the cosmic 
microwave background and the large scale structure of 
galaxies it has been possible to constrain masses of stan- 
dard model neutrinos. The bound can be derived because 
massive neutrinos contribute to the cosmological matter 
density, but they become non-relativistic so late that any 
perturbation in neutrinos up to scales around the causal 



horizon at matter-radiation equality is erased, i.e. the 
kinematics of the neutrino mass influences the growth of 
structure in the Universe. Quantitatively, neutrino free 
streaming leads to a suppression of fluctuations on small 
scales relative to large by roughly AP/P ~ — 8Q, V /Cl m 
2J. The density in neutrinos is related to the number of 
massive neutrinos and the neutrino mass by 



n v h 2 = 



Y^, m v _ N v m v 
93.2 eV ~~ 93.2 eV 



(1) 



if all neutrinos are assumed to have the same mass, h 
is the Hubble parameter in units of 100 km s _1 Mpc -1 . 
Such an effect would be clearly visible in LSS measure- 
ments, provided that the neutrino mass is sufficiently 
large, and a likelihood analysis based on the standard 
ACDM model with neutrino mass as an added parame- 
ter in general provides a bound for the sum of neutrino 
masses of roughly ^ m„ < 0.5 — 1 eV, depending on 
exactly which data is used ■ 

This should be compared to the present laboratory 
bound from 3 H beta decay found in the Mainz experi- 
ment, m„ e = (Eil^ei| 2 m?) 1/2 < 2.3 eV '% It should 
also be contrasted to the claimed signal for neutrino- 
less double beta decay in the Heidelberg-Moscow ex- 
periment 0, which would indicate a value of 0.1-0.9 
eV for the relevant combination of mass eigenstates, 
m ee = J2j ^ej m vj ■ Some papers claim that the cos- 
mological neutrino mass bound is already incompatible 
with this measurement. 

However, as with almost all likelihood analyses of pa- 
rameters beyond those in the simplest ACDM model, it 
is based on a relatively limited parameter space. Us- 
ing a much more complicated model with a non-power 
law primordial power spectrum it is possible to accomo- 
date large neutrino masses, provided that Type la super- 
nova data is discarded Q. Here, instead, we provide a 
very simple (from a cosmological point of view) exten- 
sion of the standard ACDM model in which the dark en- 
ergy component, X , is represented by a fluid with a more 
general equation of state, Px = wpx- For simplicity we 
take w to be a constant. Such models have been studied 
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extensively in the literature |8j and many of them are 
motivated by scalar field models, or models with mod- 
ified gravity on large scales. Most dark energy models 
have w > — 1, but it is possible to construct models with 
w < — 1 from non-standard kinetic terms in string theory 
or from models with modified gravity [jj. 

In fact, it is by now customary in almost all papers on 
cosmological parameter fitting to allow w < — 1. Here, 
we also allow w to take values below —1 in the so-called 
phantom energy regime, and perform a standard likeli- 
hood analysis for this model. As will be seen below, this 
relatively simple extension of the ACDM model vastly 
decreases the precision with which present observations 
can constrain neutrino masses. 

Likelihood analysis — We have performed a likeli- 
hood analysis using the most recent observational data 
from cosmology. We use the CMB measurements by the 
Wilkinson Microwave Anisotropy Probe (WMAP) satel- 
lite , the galaxy power spectrum provided by the Sloan 
Digital Sky Survey (SDSS) collaborati on id and the 
type la supernova data from Riess et al [llj . 

As parameters in the likelihood analysis we use a stan- 
dard flat, dark energy dominated model with the fol- 
lowing free parameters: SIcdm, the CDM density, fib, 
the baryon density, Hq, the Hubble parameter, n s , the 
scalar spectral index of the primordial power spectrum, 
r, the optical depth to reionization, Q, the normaliza- 
tion of the CMB spectrum, 6, the bias parameter, and 
w = Px I px ! the equation of state of the dark en- 
ergy. Finally we also use the contribution of neutrinos, 
£l v h 2 = tStt^t. From here on we assume three neutrinos 

u 93.2 ev 

with degenerate masses so that f2„/i 2 = 93 ™g V ■ From 
the flatness criterion this gives the dark energy density 
as Qx = 1 — ^cDM — ^6 — As additional data we 
use the data from the Hubble Space Telescope (HST) Key 
Project on H which give H = 72±8 kms^ 1 Mpc" 1 [llf . 
We allow the bias parameter b to vary freely. While it 
is possible to obtain stronger bounds on neutrino masses 
by adding information on the bias parameter, this infor- 
mation is most likely dominated by systematics and the 
errors correspondingly difficult to quantify. 



In Fig. 1 we show results for A% 2 



X 



Xo as a func- 



tion of neutrino mass, marginalized over the other rel- 
evant parameters (7 if w is kept constant, 8 otherwise). 
The marginalization was performed using a simulated an- 
nealing algorithm The normalization constant xl is 
taken to be that of the best fit model in each case. For 
a free w the best fit model has Xo = 1625.5 for 1513 de- 
grees of freedom (x 2 /d.o.f = 1.074), and for fixed w = — 1 
the best fit has xl = 1626.9 for 1514 de grccs of freedom 
(x 2 /d.o.f = 1.075). Note that for 1-dimensional parame- 
ter constraints we take the 95% C.L. to be A% 2 = 4. For 
2-dimensional parameter estimates we take the 68% C.L. 
at A x 2 = 2.31 and the 95% C.L. at A x 2 = 6.17. 

From the figure, and from table I (the bottom line), it 
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FIG. 1: A\ 2 as a function of neutrino mass for various dif- 
ferent data sets and parameter assumptions. The dashed line 
is for a fixed w = -1, using WMAP, SDSS, HST, and SNI-a 
data. The full line is for a free w with the same data. The 
horizontal lines show A\ 2 = 1 and 4, corresponding to 68% 
and 95% C.L. respectively. 



TABLE I: The 95% C.L. upper bound on the sum of neutrino 
masses from present cosmological observations. 



Data used 


J2 m » 


WMAP+SDSS+HST+SNI-a 


1.48 eV 


WMAP+SDSS+HST+SNI-a (fixed w) 


0.65 eV 



is clear that with w — — 1 we retrieve the known result, 
£m„ < 0.65 eV (95% C.L.), published in the litera- 
ture for the same data 0]. However, as soon as w is 
allowed to vary the mass bound degrades tremendously, 
and in fact allows very high neutrino masses. For the 
data used here the bound is m v < 1-48 eV (95% C.L.). 
Note that if neutrino masses are close to saturating this 
bound they will be easily detectable by the KATRIN ex- 
periment, which has a projected sensitivity of 0.2 eV for 
the effective electron neutrino mass 15], corresponding 
roughly to ~ 0.6 eV. It is also completely com- 

patible with the claimed detection of a non-zero neutrino 
mass by the Heidelberg-Moscow experiment, with a best 
fit around m ee ~ 0.1 - 0.9 eV @. 

In Fig. 2 we show the analysis in a grid for both ^ m v 
and us. From this it can be seen that there is an al- 
most perfect degeneracy between these two parameters, 
an increasing m v can be compensated by decreasing 
w. While for low neutrino masses a cosmological con- 
stant (w = — 1) is allowed, for high neutrino masses only 
dark energy models in the phantom regime (w < — 1) are 
allowed. 
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FIG. 2: 68% and 95% allowed contours as a function of neu- 
trino mass and dark energy equation of state using WMAP, 
SDSS, HST, and SNI-a data. 

The reason for the degeneracy is that when il„ is in- 
creased, Q m must be increased correspondingly in order 
to produce the same power spectrum. However, when 
w = — 1 an increasing fi m quickly becomes incompati- 
ble with the supernova data. This can be remedied by 
simultaneously decreasing w because of the well-known 
Q m , iv degeneracy in the supernova data. This effect can 
be seen in Fig. 3: If w is allowed to vary freely, fi m can 
take very high values without being inconsistent with the 
supernova data, because in this case the w = — 1 upper 
bound on f2 m does not apply. As soon as the model with 
fixed w = — 1 approaches am v ~ 0.7 eV and the best 
fit Vl m crosses the Riess et al. bound the model 
becomes strongly disfavoured. 

This model gives probably the simplest example of how 
to relax the very stringent cosmological neutrino mass 
bound. Other means have been discussed in the litera- 
ture, such as broken scale invariance ,7| or mass- varying 
neutrinos , but this is by far the simplest scenario yet 
discussed. 

Future constraints - Since present data clearly 
do not give very stringent constraints on ^ m v in 
the presence of phantom energy, it is worthwhile 
discussing whether future data will be able to break 
the degeneracy. For that purpose we have per- 
formed a Fisher matrix analysis based on a reference 
model with the same free parameters as in our fit to 
present data (Ocdm^ 2 , w x, ^bh 2 , fix, n s , r, Q, b, J2 m v ). 
The parameters in the reference model are 
(0.1225, -1, 0.0245, 0.7, 1, 0.05, 1, 1, 0.05 eV). Note 
that we again assume three massive neutrinos with com- 
pletely degenerate masses. While this approximation 



FIG. 3: The value of Sl m for the best fit models, as a function 
of ^m„. The curve labels are the same as in Fig. 1. The 
horizontal (red) lines are the 2a bounds from the present Riess 
et al. supernova data |Tl|l for the case of w = — 1. 

is valid as long as J2 m f 0.15 eV, it breaks down 
for smaller masses. This introduces a relatively small 
numerical error, and does not qualitatively change the 
results (see also [21| for a discussion of this point). 

The Fisher matrix analysis is based on the second 
derivatives of the likelihood function around the reference 
model and allows for an estimate of the la error bars on 
measurements of all the relevant parameters in a given 
experiment, assuming that errors are Gaussian. For the 
future CMB data we assume an experiment which mea- 
sures both temperature and iJ-type polarization. It is 
assumed to be limited by cosmic variance up to I = 2500 
for temperature and I — 1500 for polarization, i.e. an ex- 
periment slightly better than the upcoming Planck satel- 
lite 01 (hereafter we call this hypothetical experiment 
Planck+). We neglect B-type polarization even though it 
will most likely be measured by Planck. Adding informa- 
tion on i?-type polarization would not significantly alter 
the results. For a future LSS survey we assume an effec- 
tive survey volume of V — §ttA 3 with A = 1000 hr 1 Mpc 
(see [3 

for a discussion). This should be compared to 
the SDSS-BRG survey which has A ~ 620 /i" 1 Mpc. 
We also assume that the linear power spectrum can be 
inferred with only sampling error up to a wavenumber 
of k = 0.15/i/Mpc. In table II we show the la error 
bars which can be expected on w and ^ m„ for such a 
combination of data. 

From the table it is clear that if either w or ^2 m v can 
be assumed to be fixed then a very stringent constraint 
can be obtained on the other parameter. However, as 
soon as both parameters are allowed to very freely the 
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TABLE II: Estimated la errors on w and "^2 m v for the Planck 
+ SDSS data set. The errors have been calculated from a 
Fisher matrix analysis, as described in the text. 





fixed w or ra„ 


free w and m v 




0.106 eV (w fixed) 


0.288 eV 


a(w) 


0.026 (m„ fixed) 


0.069 



error bars blow up by a factor of almost 3. The estimate 
that cr(E TO ^) ~ 0.06 - 0.07 eV ;20|, is clearly too opti- 
mistic unless additional information on w is provided (see 
also |ilj for further discussion). Such information could 
for instance come from large scale weak lensing surveys 
|22| , or from measurements by the SNAP supernova sur- 
vey satellite |23j |. 

Conversely, in order to provide a stringent constraint 
on w it is necessary to obtain prior knowledge about the 
neutrino masses. This could come from experiments such 
as KATRIN which is designed to probe the effective 
electron neutrino mass to a precision of about 0.2 eV, or 
from neutrinoless double beta decay experiments which 
in theory are sensitive to the sub-0.1 eV range. 

Conclusion — We have studied the cosmological neu- 
trino mass bound in cosmological models where the dark 
energy equation of state is allowed to take on an ar- 
bitrary, but constant value. We find that this relaxes 
the present cosmological bound on neutrino masses by 
more than a factor of two, to ^2, m v < 1-48 eV at 95% 
C.L. Furthermore, even with the much more precise CMB 
and LSS data available in the future, the degeneracy per- 
sists unless additional data from weak lensing or similar 
probes can be used to break it. 

The example provided in this paper clearly illustrates 
that while cosmological bounds on particle physics pa- 
rameters are very impressive, they are also model depen- 
dent. A cosmological neutrino mass bound cannot stand 
completely alone, it should be complemented by direct 
laboratory measurements in the 0.1 eV sensitivity range. 
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